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Abstract 

We construct a probability model seemingly unrelated to 
the considered stochastic process of coagulation and fragmen- 
tation. By proving for this model the local limit theorem, we 
establish the asymptotic formula for the partition function of 
the equilibrium measure for a wide class of parameter func- 
tions of the process. This formula proves the conjecture stated 
in Q for the above class of processes. The method used goes 
back to A.Khintchine. 



1 Introduction and Summary. 

The motivation for our research came from a conjecture stated in 
p. 462, in the following setting. 
For a given integer A^, denote by 

N 

rj = {ni, . . . ,12^) : < Uk < N, ^knk = N (1.1) 

fc=i 

a partition of into Uk groups of size /c, k = 1,2, . . . ,N and by 
VLn = {q} the set of all partitions of A^. 

We will be interested in the particular probability measure /i on fi^v 
given by 



fiN{r]) = Cn \ ^, " " ^, , T] = {ni, . . . ,nN) e Qn, 
niln2l . . . riN'. 



[1.2) 



where > 0, k = 1,2,...,N and Cjy = Cn{o,i, ■ ■ ■ ,aN) is the 
partition function of the distribution fijy : 

N=l,2,... (1.3) 

The probability measure /ia? is the equilibrium state of a class of 
reversible coagulation-fragmentation processes (CFP's)(see for refer- 
ences Isli). 
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CFP's trace their history from Smoluchowski(1916) and they have 
been intensively studied since this date. The process models the 
stochastic evolution in time of a population of particles distributed 
into groups that coagulate and fragment at different rates. The 
model arises in different contexts of application:polymer kinetics, as- 
trophysics, aerosols, biological phenomena, such as animal group- 
ing, blood cell aggregation, etc. Observe that particular choices of 
Ofc, k = 1,2, ... ,N in ( |1.2| ) lead to a variety of known stochastic 
models. For example, when = f , k = 1,2, . . . , N, (3 > 0, ( [1.2|) 
becomes the widely known Ewens sampling formula that arises in 
population genetics. 

Following PI, we view CFP as a continuous-time Markov process 
on the state space Qn- Formally, a CFP is given by the rates 
and of the two possible transitions: coagulation and fragmentation 
respectively. Namely, ipii,]), 2<i+j<Nis the rate of merg- 
ing of two groups of sizes i and j into one group of size i + j, and 
2 < i + j < N is the rate of splitting of a group of size i+j 
into two groups of sizes i and j. We consider the class of CFP's for 
which the ratio of the transition rates has the form 



t,3 -.2 <t+j < N, 



:i.4) 



where > 0,k = 1, . . . ,N are given parameters of the process. 
Owing to (|1.4|) , the condition of detailed balance holds, and, conse- 
quently, the CFP considered is reversible with respect to the invariant 
measure (|1.2|). 

Letting — oo, we will be concerned with the relationship between 

two infinite sequences {an}'^ and {c„}[f ,co = 1. 

It was conjectured in p|* that the existence of the limit 



lim > (1.5) 

a„+i 
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implies the existence of the hmit 

hm ^ > 0. (1.6) 

Apart from the fact that the conjecture is a challenging mathematical 
problem, one can see from that it also has a direct significance 
for the stochastic model in question. First, if the limit exists, 
then a variety of functional of the process( e.g., the expected values 
and variances of finite group sizes), as oo, can be explicitly 

expressed via this limit. Next, by formula (4.16) in ^ we have that 

cov{nk,ni) = a[k)a[l) — -— , 

V Cn [cnY ) 

fc^/ = l,2,...,iV, k + l<N. (1.7) 

Thus, the validity of the conjecture will imply that at the steady 
state the random variables rifc, ri;, k ^ I become uncorrelated, as 
N oo. This fact incorporates into the assumption of independence 
of sites in mean-field models, as iV — > oo, that is commonly accepted 
in statistical physics. 

Another motivation for our study is provided by a quite different field, 
known as random combinatorial structures(RCS's). The connection 
of CFP's to this field is based on the following observation made in 
Let Zi,i = 1, ... ,A^ be independent Poisson random variables 
with respective means > 0, i = 1, . . . , N. Then it is easy to see 
that the distribution /i^v admits the following representation 

TV 
i=l 

r] = (ni, . . . ,nAr) G Qn- (1-8) 



It turns out that (|1.8|) is the general form of distributions arising in 
a variety of RCS's. This is explained in [|T|, and |TD|. (Theorem 1, 
p. 96 in 0] gives a rigorous proof of this fact). The simplest example 
of a RCS is a random choice from A^! permutations of A^ objects. 
Cauchy's formula for the number of permutations having cycles 
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of length k, k = 1, . . . , N, where Y,k=i = N, tells us that the 
probability of picking a permutation with this property is given by 
( [OD with ak = k-\ k = 1,. . . ,N and Cn = I. 

Added in proofs:The conjecture was recently proved by J. Ball and 
S.Burris in" Asymptotics for Logical Limit Laws", 2001, Preprint. 
In view of this, one can translate the preceding reasoning in the 
context of the Conjecture, into the language of RCS's. 
Our paper is devoted exclusively to the study of the asymptotic be- 
haviour, as n — > oo, of the quantity c„, defined by ( |1.3| ). 
The asymptotic formula for Cn established in our paper proves the 
Conjecture for a wide class of parameter functions a : a{k) = a^, k = 
1,2,.... We mention also two other applications of our result related 
to global characteristics of CFP's(=RCS's). 

(i) For a given n, denote by Vn the mean value of the total number of 
different groups at the equilibrium of CFP (=components in a RCS). 
It follows from (4.15) in H that 

" Cn-k 

Vn = Y.^k-^- (1-9) 

(ii) Denote by poo the probability at the steady state of the creation 
of a cluster (= component in a RCS) of infinite size. It was shown in 
[|], p. 462, that the condition 

lim {Vn - Vlanl^i) > 0, (I-IO) 
n^oo L J ' 

for some < a < 1 is sufficient for > 0. (Here [•] is the integer 
part of a number.) 

Thus, with the help of the asymptotic formula for c„, one can reveal 
the asymptotic behaviour of as n — oo and, consequently, find 
the limit in ( |1.10| ). The latter will answer a question which is common 



in statistical physics. 

Also we want to point out that determining the asymptotic properties 
of partition functions for interacting particle systems is a difficult 
mathematical problem widely discussed in statistical physics (see e.g. 
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2 Description of the method and 
a sketch of its history 

We assume a„ > 0, n = 1, 2, . . . and that the following limit exists: 

lim^^:=i?, 0<i?<oo. (2.11) 

a„+i 

Thus, the power series in 

oo 

S{x) := ^a„x", (2.12) 

n=l 

has radius of convergence R and it converges in the complex domain 
D C {\x\ < R}. 

Then (see |^), g(x) = e^^^\ x G D is the generating function for the 
sequence {cn}^^ defined by (|1.3|). Namely, 

oo 

^(x) = e^(^) = 5: c„a;", x E D, (2.13) 

n=0 

and moreover, the series (|2.12|) and (|2.13|) converge in the same do- 
main D. 

The method we use here for deriving the asymptotic formula for Cn 



goes back to A. Khintchine's pioneering monograph |TT[. In ||TT 
Khintchine developed the idea of expressing of values of quantum 
statistics via the probability function of a sum of correspondingly con- 
structed independent integer- valued random variables. Subsequent 
implementation of the local limit theorem resulted in the method of 
the derivation of asymptotic distributions of quantum statistics. In 



1^ this method was systematically applied to systems of photons 
and some other models. The method was further developed by A. 
Postnikov and G. Freiman, (see for references |]I1|) who applied it to 
analytic number theory. In particular, G. Freiman formulated a local 
limit theorem for some asymptotic problems related to partitions. A 
general scheme for the derivation of asymptotic formulae for these 
kind of problems was outlined by G. Freiman and J. Pitman in 
( for references see also H.) 
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A similar approach, also based on the implementation of the local 
limit theorem, has been independently developed for the last fifteen 
years in the theory of RCS's. A very good exposition of this direc- 



tion of research is given in the recent monograph |jTO[ by V. Kolchin. 
We will explain briefly the basic difference between the problem ad- 
dressed in the present paper and those in |T^. In the context of the 



generalized scheme of allocation that encompasses a variety of RCS's, 
5* and g are the generating functions for, respectively, the total num- 
ber of combinatorial objects of size n and for the number of such ob- 
jects possessing a definite property. In this setting it is assumed that 
the expression for the function g is known explicitly. Based on this, 
the combinatorial quantity in question is expressed via the probabil- 
ity function of a sum of i.i.d. discrete random variables, distributed 
according to a probability law that depends on the given values of 
{c„}[f . Such scheme is applicable for example, for investigation of 
the asymptotic of the number F^^n of all forests of nonrooted 
trees having n vertices, in which case c„ = (n!)~^n"~^, n = 1,2,.... 
First, w.l.o.g. we assume throughout the paper that the common 
radius of convergence of the series ( ^.121 ) and ( f^.l3D equals 1. This 



is due to the fact that taking in ( p..3|) Oj = Waj, j = 1, . . . , A^, it 
follows from (|1 . 1|) that cn = cn- The above assumption makes 



the lim„^oo if it exists, equal 1. 

Our starting point is the following representation of c„. 

Lemma 1 

Cn = e"'^/^ niE^^^^^ Jxe-2™"^ia, n=l,2,... 

(2.14) 

for any real a. 



Proof: 

It follows from ( |1.3D that c„ depends only on Oi, . . . , a„, which means 
that the first n -|- 1 terms of the Taylor series expansions of the two 
functions 
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= e^^""^ and 5f„(x) := eS"=i x G -D 



(2.15) 



are the same, i.e. Ck = Ck,n, k = 0, . . . ,n, where {ck,n}kLo is the 
sequence related to the function Qn- For a fixed n, the series expansion 
of the function gn{x) converges for all x. So, we can set 



X = e-'^+2-*", (2.16) 



for some real a and a. 
Then we have 



(2.17) 



The last equality is due to the fact that 



1 . 1, if m = , , 

e2™"^rfa=<! ' 2.18 
0, if m ^ 0, m G ^. 



Finally, substituting 



1=1 1=1 \k=0 J 



and ( gA6D in the LHS of (^TzD we get the claim. 



Our next step will be to give a probabilistic meaning to the expression 
(ra) for c„. 



We introduce the following notations. 

n 

Sn{x) = Y.aix\ (2.20) 



k\ exp (a/e" 



Pik = TTZZTTTZZ^^ / = l,...,n, A; = 0, 1, . . . (2.21) 



oo 



y.^(«) = EP^A.e2™'^ ^^^^ (2.22) 

fc=0 
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^ia) = l[ipiia), aeR. (2.23) 

1=1 

Now ( p.l4| ) can be rewritten as 

c„ = e'^'^e^"^'"'^ (^(a)e-2™"cia. (2.24) 
Jo 

The fact that for a given /(!</< n), pik, = 0, 1, . . . is a Poisson 
probabihty function with parameter aie~'^\ suggests the following 
probabilistic interpretation of the integral in the RHS of (|2.24| ). 
Let Xi, . . . ,Xn be independent integer- valued random variables de- 
fined by 



Pr{Xi = lk)=m, l = l,...,n, A; = 0,1,... (2.25) 

Then <f{a) defined above is the characteristic function of the sum 
Y = Xi + . . . + Xn and we have 

r (^(a)e-2"^'""rfa = Pr(Y = n). (2.26) 
Jo 



Now ( |2.24| ) can be viewed as an analog of the aforementioned Khint- 
chine's representation for c„. 

It is well- known from the classical theory of limit distributions of 
sums of independent integer- valued random variables that, under cer- 
tain conditions on distributions of the variables, a local limit theorem 
is valid. 

In our subsequent study, the free parameter a will be taken depending 



on n : a = an- By (|2.25|) and (|2.21| ) this means that the probability 



law of each of the n random variables Xi, I = 1,2, . . . ,n depends 
on n. Therefore, to compare with the classical case, we will be deal- 
ing here with a triangular array of random variables. For this case, 
general necessary and sufficient conditions for validity of the local 
limit theorem are not known. For some cases results in this direc- 



tion were obtained in [13| and H]. In the first of these two papers a 



sufficient condition was established (see (|]T3[, Theorem 2, condition 
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Ill) in the case of an array of general lattice random variables. It can 
be verified that this condition ( which can be viewed as a version of 
the celebrated condition of asymptotic uniformity) fails for the class 
of parameter functions a considered in our paper. The second pa- 
per studies exclusively the case of the triangular array of trinomial 
random variables. 

Most of this paper is devoted to the proof of the local limit theorem 
in the above setting. 

Namely, we will demonstrate that under certain conditions on the 
parameter function a 

Pr(F = n) ~ (27r5^)-^/2g-(M„-n)V2B^^ 72 ^oo, (2.27) 

where M„ = EY and Bl = VarY. 



3 Proof of the local limit theorem 



In order to prove ( p.27| ) we have to find the asymptotic formula, as 
n —>■ oo, for the integral in the RHS of ( |2.24| ). We will denote in the 
sequel, by 7, 7^, 2 = 1,2,... constants. 

First, we obtain the explicit expressions for the quantities M„ and 
Bl 



(|2.25|) and ( |2.21| ) say that / ^Xi, 1 = 1,..., are Poisson(a;e '^') ran- 



dom variables. So, we have 

EXi = laie-"\ l = l,...,n, (3.28) 

VarXi = l^aie-^"", l = l,...,n. (3.29) 



This gives 



Mn = Y.^aie-'% n = l,2,... (3.30) 



1=1 

n 



Bl = ^/%e"'", n = l,2,... (3.31) 

1=1 
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It follows from the preceding discussion that the representation ( [2.24 ) 



holds for any real a. Our next result shows that a can be chosen so 
that the exponential factor in the RHS of (|2.27| ) equals 1, for any 
n = 1,2, 



- 5 -"5 



Lemma 2 

The equation 

n 

Y^laie~^" = n (3.32) 
1=1 

has a unique solution a = an, for any n = 1,2,.... 

Proof: The assertion follows immediately from the assumption ai > 
0, / = 1,2,... ■ 

Remark The above choice of the free parameter a makes the prob- 
ability of the event {Y = n} large, as n —>■ oo. The same idea is 
widely used for approximation of RCS's by independent processes 
(see for references |1|). In statistical physics such a way of choosing 
a free parameter for estimating averages is known as Darwin-Fawler 
method developed in the 1930's (for references see, e. g. 0) 
It follows from ( p.32| ) that if the series lo-i converges, then cr„ < 
for sufficiently large n, while in the opposite case the sign of 
depends on the behaviour of 5*^(1), as n oo. However, in both 
cases the following basic property of a^, n = 1, 2, . . . holds. 

Lemma 3 

Let 

lim ^ = 1. (3.33) 



n— »oo 



Then lim^^oo = 0. 
Proof: 

By the definition of (J„, 



n-^Y.^aie-^"" = 1, n=l,2,... (3.34) 

1=1 

Denote bi = lai > 0, I = 1,2 ... . Based on ( p. 331 ) 
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let N = N{e), e > be s.t. 



1 - e < ^ < 1 + e, for all / > A^. (3.35) 
k 

Consequently, 

(1 - ey-^bN <bi<{l + ey-^bN, I > N. (3.36) 

Now suppose that limfc_,oo o"nfe = o", for some subsequence rik — > +oo, 
a.s k oo, where |cr| < oo. 

Let first — oo < a < 0, then taking < e < 1 — e°"/^ we have 

1=1 l=N 

as /c — >• OO, (3.37) 

since 

(1 - e)e"'^"fc > e^~'^"fc -> e'"^'^ > 1, as k ^ oo. 
If now < a < oo, then for < e < e'^/^ - 1, 



1=1 

\Z=1 l=N 

as k ^ oo, (3.38) 



:i K^5^e-'-.. + (l + e)-^6^^((l + e)e--^) ^0 



since m this case e''^/^ < ^ Bq^^ (|337|) and ( g^SD contradict (|3]3|) 



which implies that limsup^^oo^™ = hminf„^ooC"n = 0. ■ 

In what follows we will assume that cr = cr„, as defined by ( |3.32| ). Our 
next lemma provides the expression for the integrand in the LHS of 
( |2.26|) for small values of a. 



Lemma 4 

For a fixed n and a = an, 



¥p(a)e-^™" = exp i^-2TT^a^Bl + 0{a^ps)) , as a ^ 0, (3.39) 
where ps = psi^n) = Er=i l^aic'^"", n = 1, 2, . . . . 
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Proof: By - CT) , 



k 



aeR (3.40) 

and 

<f{a) = exp a^e-'"" (e^™' - l)^ , a G i?. (3.41) 



Finally, substituting in (p.41|) the Taylor expansion (in a 



e^™' -1 = 27rm/-27rV/2 + 0(a¥), as a ^ 0, (3.42) 



that holds uniformly for / > 1 , and making use of the definition ( p.32|) 
of (J„, proves ( p.39|) . 

Observe that uniformity (=the constant implied by the term 0{a^P) 
in ( |3.42| ) does not depend on I) is due to the fact that for any real a, 



da' 



e'™'| < (27r)T, s = l,2,.... 



Now we are prepared to deal with the central objective stated in the 
beginning of this section. Denote 

T = T(n) = /V(a)e"^™"cia, n = l,2... (3.43) 
Jo 



The integrand in ( p.43| ) is periodic with period 1. So for any < 



tto ^ 1/2, the integral T can be written as 



T = Ti + T2, (3.44) 

where Ti = Ti(q;o; n),T2 = T2{ao; n) are integrals of the integrand in 
(|3.43| ) over the sets [— ao, +cto] and [—1/2, — Q;o]U[a;o, 1/2] correspond- 
ingly. Following the idea of ||^, [§], we will first show that for an ap- 
propriate choice of ao = aoi^^) the main contribution to T, as n ^ 00 
comes from Ti. Then, estimating Ti, under ao = oiQ{n), n — 00 we 
will get the desired asymptotic formula ( p.27| ). 
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It is clear from Lemma ^ that the asymptotic behaviour as n ^ oo 
of the integral Ti is determined by the asymptotics of the three key 
parameters cr„, and ps. 

First we address the problem for the class of parameter functions a 
of the form 

ttj = j = 1, 2, . . . for a given p > 0. (3.45) 

It follows from the definition of the above three parameters that in 
the case considered the problem reduces to the estimation of sums of 
the form Z)j'=i j^e"'^-', k > —1, as n — > oo. 
To do this we apply the integral test for the function 

f{u) = (aM)'=e-"", M > 0, a > 0, A; > -1. (3.46) 

In the case k > the function / is strictly increasing on [0, ka~^] 
and is strictly decreasing on [ka~^, +oo). So, applying the integral 
test separately on each of the above intervals we have in the case 
considered 

/ f{u)du+jJ{ka-')<Y.fij)< f{u)du+j2f{ka-'), (3.47) 

where the constants 71,72 depend on k only and f{ka~^) = k^e~^. 
If now — 1 < /c < 0, then the function / is strictly decreasing on 
[0; +00) and the integral test gives 

/ f{u)du<Y.fU)< / f{n)du. (3.48) 
Ji ~{ Jo 

Note also that for any a,b > 

rn rn rna 

/ f[u)du = a'^ / u'^e-'^'^du = a"^ / z''e-'dz. (3.49) 

Jb Jb Jba 



Then, combining ( p.46|) -( p.49|) we arrive at the desired asymptotic 
estimate: 



n 



E^'^'e"'-^, as 0<(T^0, n-^oo, 
j=i ^ 
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and liminf no" > 0. (3.50) 

n— ►oo 

In particular, if ^ oo, then the constant 7 in ( ^.501 ) can be found 
exphcitly: 

tf^-" ~ (3-51) 
where F is the gamma function. 

Further we will write •{n) x if there exist positive constants 71, 72, 

s.t. 7in" < •{n) < 72^", for all sufficiently large n. 

Extending (|3.45| ) we consider now the class of functions a satisfying 



^, as j ^ 00 for a given p > 0. (3.52) 



An obvious variation of the preceding argument gives in this case the 
following analog of ( |3.50|) : 



" . 1 

J2 X — y, as 0<(T->0, 



,=1 



and liminfncT>0, / > 0. (3.53) 

n^oo 

This immediately implies 
Lemma 5 

Let the function a obey ( \3.53i ). Then, as n ^ 00, 



an^n p+i (3.54) 
Bl X niTT (3.55) 

and 

P+3 

P3 X np+i . (3.56) 

Proof: 

By the definition of cr„, Lemma 3 and ( p.53| ), in this case 

n X , as n ^ 00 (3.57) 



p+i' 

On 



and, consequently. 
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(T„ X n p+i, as n — >• oo. (3.5J 



Now the last two assertions follow from ( p.53|) . 



At this point we are prepared to estimate the integral 

Ti = Ti{ao;n) = (/?(a)e-2™"da. (3.59) 

■J —an 



Lemma 6 

P+2 

Let the function a obey ( \3. 5^ ) and ao = an'^ log n, p > 0. 



Then 

Ti{ao;n) ^ {27TBly^/^, as n ^ oo. (3.60) 

Proof: 

By ( |3.54| ), ao — > 0, as n ^ oo. So, making use of ( ^.39| ) we obtain 

f-ao 



/CtO / s 

exp (—2TT'^a^B'^ + O^a^ps)] da, as n ^ oo. (3.61^ 
-an ^ ' 



Further, under the condition ( p.52| ), it follows from Lemma 5, that 
the above choice of ao provides the following two basic relationships 
between the parameters B\ and ps : 



lim an-B„ = lim log n = +oo (3.62) 

and 

£ 

lim anPs = lim an log^n = 0. (3.63) 

n^oo " n^oo 



Since lim„^oo a^Ps = for all a G [— ao, ao], ( p.61|) and (|3.62| ) imply 



Ti ~ /" exp [— 27r^a^-B^] lia = 

J-«o 

]^ r2TTaoB„ ^ \ 

— - / exp(-— )(iz , as n ^ oo. (3.64) 



Taking ao as in Lemma 6, we write 

T2 = Ta,! + T2,2, (3.65) 
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where 

T2,i = r2,i («o; n) := j^'^ ¥;(«)e-2™"rf« (3.66) 

and T2^2 is the integral of the same integrand, but over the set 
[-1/2, ao]. 

In view of ( p.65| ) and the fact that a) = V^(a), a G -R, the rest of 
this section is devoted to estimation of the integral T2^i, as n oo. 
Our starting argument will be the same as in It follows from 
that 



\ip{a) I = exp aje~^"^ sin^ vraj j , a e R. (3.67) 

Denote by [x] and {x} respectively the integer and fractional parts of 
a real number x and the distance from x to the nearest integer, 
so that 



, „ , {x}, if {x} < 1/2 ^ ^ 

^ l-{x}, if {x}>l/2. 
We will make use of the inequality 

sine>-^, O<0<-. (3.69) 

TT 2 

Since sin^ vrx = sin^ 7r||x|| for any real x, it follows from (|3.69| ), ( p.68|) 
that for all real x 



sinVx>4||xf. (3.70) 
Hence, in view of ( ^.67] ) and ( 3.66| ) we have to estimate the sum 



V^{a) := 5Zaje"^'""||ajf , ao < a < 1/2. (3.71) 

Lemma 7 

Let the function a obey (\3.5^ ). Then 



\4(a) > 7 log^ ra, a G [ao, 1/2], as n oo, (3.72) 
where 7 > 0. 
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Proof: Take ai = ai{n) = On > «o(^)) as n ^ cx) and split the 
interval [ao! 1/2] into two disjoint subintervals l\ := [oq, ai] and I2 '■= 
(ai, 1/2]. We plan to prove the assertion separately for a G /i and 
a G /2. 
Interval Ii. 

It is clear from the definition (|3.68|) that 



= aj, j < 77—^ " ^ ^1- (3-73) 

In view of the fact that, by Lemma 3, {2ai)^^ 00, as n 00, 
while (q;i)~V„ = 1, n = 1,2, we apply ( p.53|) with / = 2 to 
obtain 



Vn{a)>al J2 a,jV^"" x^=log2n, a G Ji. (3.74) 

l<i<(2oi)-l 

Interval J2. 

For a given integer n and a given a G /2 define the set of integers 

Q{a) = Q{a;n) = {l<j<n: \\aj\\ > 1/4}. 
It is clear that 

Q(a) = {1 < j < n : k + l/A < aj < k + 3/4, A; = 0, 1, . . . , } = 



UU 'Qk{a), (3.75) 

where (^^(a) denotes the set of integers {j : ^|il l£ j ^ ^i^}- 
Observe that for any a G /2 and k > the set Qfc(a) is not empty, 
since in this case — > 1. 

This yields the following estimate of the sum V^(a), a G /2 : 



Vn{a) > 1/16 J2 cLje'^"" = 1/16 J2 J2 a^e"^'"". (3.76) 

jeQia) k=0 jeQfe(a) 

We now assume that the asymptotic inequality (|3.52|) holds for all 
j > N. This means that ( p.52| ) is valid for all j G Qfc(«) whenever 



k > max{0; ^^^} =-. K{a). 
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We agree, with an obvious abuse of notation, that for a real u, Quicn) 
is the interval 

Au + 1 4m + 3 

^ 4a ' 4a ^' 
Observe that for all sufficiently large n, we have < a~^cr„ < 1, 
a G I2, while nan 00. 

Applying now the integral test to the double sum in the RHS of ( p.76|) 
gives 



Vn{a) > 7 / du vP-^e-'"'"dv = 

JK{a) JveQu{a) 

-1 4an — 3 

7^ / ' du [ yf-^e-'^dv > 

an JK{a) Jve{cr„Qu{a)) 

7^ / ° v^~^e~'"dv I" du^ 

an JKi{a) an 

p > 0, as n — > 00, (3.77) 
where we denoted Ki{a) = ^^^^^^Cn and 

4m + 1 4m + 3, 



CTniQuia)) = [ar, 



4a ' 4a 



Note that the last inequality in ( ^.771 ) is obtained via the change of 
the order of integration. Finally, (|3.77|) , (|3.54|) and (|3.74| ) prove the 
claim. ■ 

The last statement of this section is the desired local limit theorem. 



Theorem 1 

Let the function a obey 52[ ). Then 



T = PriY = n) ~ , ^ , as n ^ 00. (3.78) 
l2TiBl 



Proof: 

By Lemma 7 and ( p.66| ), ( p. 671) , ( p.65|) , we have 



< g-27iog2„^ asn^cx), (3.79) 
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where 7 > 0. In view of ( |3.64| ), and ( |3.44| ) this proves ( |3.78| ). 



We provide now an extension of the field of vahdity of the above local 
limit theorem. 

We agree to write n^^ ^ •{n) ^ n^^, < f3i < f32, if there exist 
positive constants 71,72, s.t. 7in^i < .(^^^ < ^gn^^, for all sufficiently 
large n. 

For given < pi < p2 define the set J^{pi,P2) of parameter functions 
a = a{j), j G -R"*", obeying ( |3.33| ) and the condition 



f'~' ^ a{j) ^ f'-\ 0<pi< P2. (3.80) 

Corollary 1 

For an arbitrary p > and < e < | the local limit theorem ( p. 781) 

is valid for all parameter functions a G + e;p). 

Proof: 

It is clear from the preceding results that for all a G J^{pi',P2), we 
must have, as n — 00, 



n w+i ^ o-„ ^ n w+i (3.81) 

Therefore, setting, as in Lemma 6, ao = (-Bn)^^ log n, gives 

«^P3<73(log^n)a„^a-(-+^) = 
73(log n)an ^ — 0, as n — > 00, (3.82) 
provided 3pi — 2p2 > 0. 

Thus, in the case a G J^{pi;P2), where pi,P2 : ^ > P2 > Pi > 0, 
Lemma 6 is valid. The proof of Lemma 7 for this case goes along the 
same lines, with an obvious replacement of (|3.77|) by 



4r ^ ^-(«) ^ « e ^2. (3.83) 
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Combining these results proves the vahdity of the local limit theorem 
for the class of functions a in our statement. ■ 
For our subsequent study we will need the following extension of 



Corollary 2 

Under the conditions of Corollary 1, 

PriY = n + h) ^ Pr(Y = n) ^ -^=, as n oo, (3.84) 

for a fixed real h. 
Proof: By 



rl/2 

Tih- n) := Pr{Y = n + h) = / (/?(a)e~2™("+'')rfa, (3.85) 

J~l/2 



where the characteristic function (p{a) is given by ( p. 221 ) and ( |2.23| ) 
By Lemma 4, we get 



y?(a)e-2™("+^) = exp [-27r^a^Bl - 2mah + ©(a^s)) , 
as a — 0, (3.86) 

where ps is defined as in ( p.39| ). Next, let ao = ao('^) be as in 
Lemma 6. Denote 

(^(a)e-2™("+^)da. (3.87) 
Lemma 6 and (|3.86|) imply for a fixed h & R 

Ti{h; n) =, as n — > cx). (3.88) 

127131 



Now it is left to estimate the integral 



f2{h-n) ■= r^V(a)e-2™("+'*)da + / "V(«)e"'™^"+''^c^«. 

Jao J-1/2 

(3.89) 

Since the function ^{a) here is the same as in ( p. 67 ), the estimate 



( |3.79|) is valid also for T2(/i; n), which together with (|3.88| ) proves the 



statement. 
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4 The asymptotic formula for 



By virtue of ( 2.24|) and Corollary 1, we obtain the following asymp- 
totic formula for c„ valid for all parameter functions a G -^(^ + 
e;p), p>0, e>0: 



: exp nan + X! '^i^' 



as n — >• cxD. 



(4.90) 



Example. Let aj ^ ^, p > 0, as j — > cxd,. Then by Lemma 5, 
( |333|) and ([CTI ), 



log c„ X 2np+i — -log27r — -^og n, as n — > oo. (4.91) 

In particular, if Oj = p > 0, j = 1, 2, . . . , then using ( |3.51| ), 



the constants in (|3.54|) , (p.55|) can be found explicitly and we obtain, 

as n oo, 

/ n 



1 

' p+i 



On ~ 



r(p + i) 



Bt 



n 



V+2 
P+1 



r(p + i: 



r(p + 2) 



(4.92) 
(4.93) 



nan ~ (r(p + l))f+i 



(4.94) 



j^a.e-^'^- r^V{p)[-^^^^y^' =p-i(r(p+l))^n^. (4.95) 



Hence, ( [4.90| ) gives, as n ^ oo 



logc„ ~ A(p)np+i - ^log27r - ^ ^ ^ (log n - log T{p + 1)) - 
2 2p + 2 

l/21og r(p + 2), (4.96) 

where 

A{j)) = {l+p-^){V{j)+l))^^ . 
Remark. For the two cases a,- = const and a,- = j, j = 1, 2, . . . the 



first (= the principal ) term in the asymptotic formula (f4.96|) was 
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obtained in by solving for large n the corresponding difference 
equations (|4.101|) below. Note that this approach is not applicable 



even for the class of parameter functions aj = ^, p > . 

With the help of ( [4.90|) we are able to address the question on the 
validity of the Conjecture stated in the Section 1 (see (pT5| ) , ([L^ )). 



Assertion. The conjecture is valid for all parameter functions 

aGJ^(f + e;p), p>0, e > 0. 

Proof: 

It is clear from ( |1.3|) that the expression for c„+i can be written in 
the following way 

Cn+l = + Cn+1, (4.97) 

where c„+i can be viewed as the value of c„+i, (given by (p..3|)) when 
a„+i = 0. Next recall that the representation ( |2.24| ) is true for any 
real a and all = 1, 2, . . . . We now apply this representation for the 
particular case ctn+i = 0. By (|2.21|) and ( |2.25|) , in this case Pr{Xn+i = 



0) = 1. So, combining (|2.24|), (|2.26|) and taking a = Un gives 



Cn+i = exp \{n + + a^e"^''" j Pr{Y = n + l), (4.98) 

where the random variable Y is defined as in Section 2. 
Substituting the expression ( [4. 981 ) in ( |4.97| ) and applying ( [4. 90] ) and 
Corollary 2, we obtain 

~ e " H , as n ^ cxD. (4.99) 

In view of the assumption ( p. 33 ) and Lemma 3, to complete the 
proof we have to show that 



lim — = 0. (4.100) 

n— »oo n 



To get this result we will implement the difference equation (see for 
references p. 460) derived from ( [2.13| ): 
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Co = 1, Ci = ai, 



(n + l)c„+i = + l)«i+ic«-.-, n = 1, 2, 

j=0 



(4.101) 



This gives for a fixed > 1, 



(j + 1) 



n-j 



O-n+l 



j=n—k 



[n + l)an+i 



(4.102) 



Consequently, by (OD and (O^) 



liminf > cq + Ci + . . . + Cfc > 1 + ai + . . . + a^, (4.103) 
for any fixed k. 

The desired conclusion now follows from the fact that if 
aeT(^ + t]p), p>0, e > 0, then 

oo 

Y,aj = oo. (4.104) 



Concluding remark 

Our asymptotic formula ( ^4.90| ) is restricted to the case 
a E ^(y + e;p)) P > 0, e > 0. The reason we require p > 
comes from the fact that the local limit theorem we proved assumes 
convergence to the normal law only. This type of convergence is 
guaranteed if ( p.62| ),( |333D hold, or, equivalently, if the parameters 
P3 and Bn obey the condition: 

lim = 0. (4.105) 

n 



Based on the reasoning preceding ( ^.281) , it is not difficult to show 



that the quantity ps has the following meaning: 



p, = Y,E{Xi-EXif. (4.106) 
1=1 
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In effect, denoting by Ui the Poisson(A/ := a/e °'') random variable, 
we have 

E{Xi - EXif = l'E{Ui - EUif = l^Xi, (4.107) 

where the last step is due to the known property of the third central 
moment of Poisson distribution (see, e.g. [|l^, p. 33). 

106|) explains that ( |4.105|) is Lyapunov's sufficient condition for 
the convergence to the normal law in the central limit theorem. To 
demonstrate that for p < the condition ( [4.105 ) fails, we consider 
the case aj = j^^, j = 1, 2, . . .. It is easy to see from ( |3.32|) that 
for this case cr„ = 0, n = 1, 2 . . . Consequently, we have, as n — > cxd, 
B"^ ~ n"^ and ps ~ n^, which gives lim„^oo = 1- 
The above discussion suggests that for p < one can expect conver- 
gence to other stable laws. The study of this case will require a quite 
different estimation technique. 
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